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Abstract We consider the "moment vanishing problem" for a general class of 
piecewise-analytic functions which satisfy on each continuity interval a linear ODE 
with polynomial coefficients. This problem, which essentially asks how many zero 
first moments can such a (nonzero) function have, turns out to be related to several 
difficult questions in analytic theory of ODEs (Poincare's Center-Focus problem) 
as well as in Approximation Theory and Signal Processing ("Algebraic Sampling"). 
While the solution space of any particular ODE admits such a bound, it will in the 
most general situation depend on the coefficients of this ODE. We believe that a 
good understanding of this dependence may provide a clue for attacking the prob- 
lems mentioned above. 

In this paper we undertake an approach to the moment vanishing problem which 
utilizes the fact that the moment sequences under consideration satisfy a recurrence 
relation of fixed length, whose coefficients are polynomials in the index. For any 
given operator, we prove a general bound for its moment vanishing index. We also 
provide uniform bounds for several operator families. 



1 Introduction 

Let / : [fl,^?] — > R be a bounded piecewise-continuous function with points of dis- 
continuity (of the first kind) 



fl = ^0<<^l <---<^p<<^p 
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satisfying on each continuity interval [^j, ^,+1] a linear homogeneous ODE 

S)/ = 0, (1) 
where J) is a linear differential operator of order n with polynomial coefficients: 

T)=pn{x)d" + --- + pi{x)d+po{x)l, d ^ deg pj <dj. (2) 

We say that such /belongs to the class The union of all such 

is the class of piecewise D-finite functions, which was first studied in ||3|. 

Any f ^ has finite moments of aU orders; 

nb 

mk{f)= / A-VWdx, ^ = 0,1,2,... (3) 
We consider the following questions. 

Problem 1. Given D and p, determine the moment vanishing index of 0^&{'£),p), 
defined as 

(7(D,p)= sup {fe: mo(/) = ---=m,(/)=0} + l. 

/G#5>(S,p),/^0 

In TheoremObelow we shall prove that the moment vanishing index is always finite. 
Consequently, the following problem becomes meaningful. 

Problem 2. Find natural families (Z 3^2! which admit a uniform bound on the 
moment vanishing index, i.e. for which 

cr(^)= sup (7 (£>,/?)< +00. 

Our main results, presented in Section @] provide a general bound for a{'£),p) in 
terms of S). As a result, several examples of families admitting uniform bound as 
above are given. The main technical tool is the recurrence relation satisfied by the 
moment sequence, established previously in ||3|. 

Our main application is the problem of reconstructing functions / G from a 
finite number of their moments. Inverse moment problems appear in some areas 
of mathematical physics, for instance heat conduction and inverse potential theory 
|[Tl|9l, as well as in statistics. One particular reconstruction technique, introduced 
in 131 and further extended to two-dimensional setting in [4|, can be regarded as 
a prototype for numerous "algebraic" reconstruction methods in signal processing, 
such as finite rate of innovation ITSl and piecewise Fourier inversion (l2]|6l. These 
methods, being essentially nonlinear, promise to achieve better reconstruction ac- 
curacy in some cases (as demonstrated recently in ||2]|6]), and therefore we believe 
their study to be important. In Section|2]below we show that an answer to Problem|2] 
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would in turn provide a bound on the minimal number of moments (measurements) 
required for unique reconstruction of any / G In essence, the results of this pa- 
per can be regarded as a step towards understanding the range of applicability of the 
piecewise D-finite reconstruction method to general signals in 3^ '3. See Section |2] 
for further details. 

Given a family ^ C l^Q), consider the corresponding family of moment generat- 
ing functions \lf (z) j^gjr; where// (z) = L^^o/Wi; (/)z^^^'. Obtaining information 
on the moment vanishing index is in fact an essential step towards studying the an- 
alytic properties of If, in particular a bound on its number of zeros near infinity 
(as provided by the notion of "Taylor Domination", see O), as well as conditions 
for its identical vanishing. In turn, these questions play a central role in studies of 
the Center-Focus and Smale-Pugh problems for the Abel differential equation, see 
ir7l [T4l and references therein. 

The moment vanishing problem has been previously studied in the complex set- 
ting by V.Kisunko lfT2l . He showed that a uniform bound a {^) exists for families 
^ consisting of non-singular operators, by using properties of Cauchy type inte- 
grals. In contrast, in this paper we consider the real setting only, while proving uni- 
form bounds for some singular (as well as regular) operator families. Our method 
is based on the linear recurrence relation satisfied by the moment sequence. Using 
this method, in Section |5] we provide an alternative proof of Kisunko's result, stat- 
ing that the moment generating function //(z) of some / G £P!^{D,p) satisfies a 
non-homogeneous ODE 

for a very special rational function Rf{z), which depends on 1) and on the values of 
/ at the discontinuities. 

In Section |6] we provide an interpretation of our main result in the language of 
Fuchsian theory of ODE. 
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2 Moment reconstruction 

We start by defining some preliminary notions. 

Definition 1. The Pochhammer symbol denotes the falling factorial 

(0/=/('-i)-----('-y + i), /gK, ;eN 
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and the expression is defined to be zero for / < j. 



Definition 2. Given D of the form (|2]i, the bilinear concomitant ( ifTTl p.21 1]) is the 
homogeneous biHnear form, defined for any pair of sufficiently smooth functions 
u {x) , V (x) as follows (all symbols depend on x): 

(«, v) = « {inv-d ip2v) + ■■■ + (-1)"-' d"-' {p„v)} (4) 

+ U' { P2V - a (p3 V) + • • • + (- 1 d"-^ [pnV) } 
+ ... 

+ m("-"-(p„v). 

Proposition 1 (Green's formula, lITTl ). Given D of the form let the formal 
adjoint operator be defined by 

s*{-} = L(-iy^n/'.w-}- 

7=0 

Then for any pair of sufficiently smooth functions u [x) , v {x) the following identity 
holds: 

V [x) (£> u){x)Ax- [ u (x) (D* v) [x) Ax = P^ (m, y) {b) - (m, v) (a) . (5) 

J a 

Tlieorem 1 (|[3l). Let f e {Tl,p) with D of the form (|2|. Denote the discontinu- 
ities of f by a = ^0 < < ■ ■ ■ < < i = b. Then the moments nij^ ~ f (x) dx 
satisf^ the recurrence relation 



LL«/j(-iy('+^);'«<W+A = e-t, A: = 0,1,..., (6) 

;=Oi=o 

where 

= - E {p^ (^;) } • (7) 

7=0 

Proof. Apply Green's formula (|5) to the identity 



J /(D/)(x)djc = 



for each y = 0, . . . , p and sum up. The result is 



^ For consistency of notation, the sequence is understood to be extended witii zeros for 

negative k. 
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it'fi')^- ^ - iy. (/,,') («,^,) M («/)} 

J=0-'ij 7=0 

The left-hand side of the last formula is precisely the linear combination of the 
moments given by the left-hand side of This finishes the proof. □ 

Now consider the problem of recovering / G i'^^p) C i^!^ from the moments 
{mo (/),..., niN (/) } (the operator 2) is assumed unknown in the most general set- 
ting). Based on the recurrence relation (|6), we demonstrate in ID that an exact re- 
covery is possible, provided that the number G N is sufficiently large. However, 
the question of obtaining an upper bound for turns out to be non-trivial, as we 
now demonstrate. 

Definition 3. Given and p, the moment uniqueness index T {'D,p) is defined by 

T (D , /7) = sup {k: ntj (/) = ntj (g) , ^ j ^ k} + I. 

In other words, given D and p, at least T (£>,/?) first moments of / G 3^31 {X),p) are 
necessary for unique reconstruction of /. 

Recalling boundedness of a{D,p) (see Theorem|3]below), we immediately obtain 
the following conclusion. 

Lemma 1. For any operator 35 and any p 

T(S),/7)<a(D,2/7). 

Proof. Let N = G {D,2p). Take/i ,/2 having p jump points each, satisfying 'Dfi = 
Q,T)f2 = on each continuity interval such that 

'"o(/i) = mo if 2) 
mNif,) = mNifi)- 

The function g = fi — fi has at most 2p jumps, and it satisfies Dg = on each 
continuity interval. The first moments of g are zero, therefore g = and thus 
/i = /2. Therefore T {D,p) <N. □ 

Consequently, in order to uniquely reconstruct an unknown / G =^ C it is 

sufficient to get a uniform bound a (,^) for the family Perhaps the most nat- 
ural choice for such families is when the parameters /?,«, are fixed. Un- 
fortunately, without making additional assumptions, the moment vanishing index of 
such families cannot be uniformly bounded. This can be seen from the following 
example. 
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Example 1. Let D,„ denote the Legendre differential operator 

D„,= (l-,Y2)-i--2x-^+m(m+l)I, 
^ ' ax^ ax 

and consider a — b ~ \ and p = 0. It is well-known that for each m e N, the 
regular solution of / = is ^„ - the Legendre polynomial of degree m. Since 
the first m~\ moments of ^„ are zero, we conclude that 

(7 (Dm) =OT 

and therefore (T (S ) cannot be uniformly bounded in terms of the combinatorial type 
of S only. 

Using the subsequent results, in Section|6]we shall in fact provide an explanation of 
this behaviour 



3 Generalized power sums 

Proposition 2. The sequence {£<:}, given by Theorem\I] is of the form 

^k = LL^j''(k),c,j, (8) 

7=0 f=0 

where each C(j is a homogeneous bilinear form in the two sets of variables 

{p,n{^j),p'„Mj),...,pi"-'\Wn=0, 

{f{^l)- f{^y), n^D- f'{ii)....j^"-'Kit)- f^'^-'K^i)}. 

Proof. Denote for convenience / (a^) ~ f = 0. Now consider the definition of 
{ea } given by (|7]l. Rearranging terms, we write 

Furthermore, using the fact that the functions {p,„ (jc)}"„=o ^'^ ^re continuous at 
each ^j, we have 
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P^[f..^){^t)-P^{f,.'^){i,j)= (9) 

+ ... 

+{/"'-" (^;) -/"'-" (^7) }p„(^i)^j-. 

Now using the definition (IDi, the claim is evident. □ 

The expression ^ for £1^ is nothing else but a generalized power sum. Let us recall 
several well-known facts about them (see e.g. [8] Section 2.3] or ifTSl '). 

Proposition 3. Let the sequence s/^ be of the form 

= L L (^)^ e (10) 

7=0 £=0 

Then it satisfies the following linear recurrence relation with constant coefficients 
of length n{p + 2) + I: 

(n(E-^,i)").,=o (11) 

where E is the forward shift operator in k and I is the identity operator 
Conversely, the fundamental set of solutions of the recurrence relation (lllb is 

{e,Mt\- ■ ■ ^ ■ ■ ■ , ^^uk^^^l- ■ ■ , ^,';r ' } • 

Corollary 1. The sequence s^ as above, which is not identically zero, can have at 
most n{p + 2) — 1 first consecutive zero terms sq = ■ ■ ■ = ^,,(p+2)-2 = 0- 

Proof. If io = • • • = ■Sh(p+2)-i = 0, then by the recurrence relation (ITTi we would 
have automatically i„(p+2) = i«(p+2)+i = • • • = 0. □ 



Corollary 2. Assume that the numbers {"^jj^^o C C are pairwise distinct. Let the 
sequence s^ be given by dlOl ), with a-priori unknown {fljj}- If si^ — Ofor all A; G N, 
then necessarily all the coefficients {fli j } are zero. 



4 Main results 



Let us now return to our main goal, namely, obtaining upper bounds on the moment 
vanishing index a (2), p) . 
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Definition 4. Given £) of the form denote for each j — Q,...,n 

def , 
«./ dj-J, 

and also 

def 

a = a (25) = max aj. 

j=0.....n 

Proposition 4. Lef / e ^S'{D,p). Then vanishing of the first {p + 2)n + a{D) 
moments of f (i.e. mq — ■ ■ ■ = ni(^p^2)n+a-l — 0) implies identical vanishing of the 
sequence {e^} defined by Theorem\l] 

Proof. Consider the recurrence relation (|6). Denote its left-hand side by /Z^ . Obvi- 
ously, since each jXf. is a linear combination of the moments, we have 

Mo = ■■■ =Mn(p+2)-l =0- 

Consequently, the corresponding right-hand sides also vanish, i.e. 

eo = --- = e„(p+2)-i =0. (12) 
The conclusion follows immediately from CoroUary[T] □ 

Now we establish our main result. 

Theorem 2. Let f G / ^ with discontinuity points 

a = ^0 < < ■■■ < < ^p+i = b. 

Assume that p,-, {^j) ^ Ofor at least one t,j as above. Then at most 

{p + 2)n + a{X))-l 

first moments of f can vanish (i.e. otq = • • • = m(p+2)n+a-2 ~ 

Proof. Assume by contradiction that the first (p + 2)n + a moments of / vanish, 
i.e. 

mo = • • • = m„(p+2)+a-i = 0. 
By Proposition|4]and Corollary |2] we immediately conclude that 

Qj = 0, j = Q,...,p + l,e = Q,...,n-l, 

where {qj} are described by Proposition |2] Now we take the concrete j for which 
Pn i^j) 7^ 0- This means that the operator J) is regular at ^j, and consequently each 
solution to £• / = in the neighborhood of is uniquely determined by the initial 
values /(^y) , . . . i^j). We claim that 
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(13) 

In this case, we would immediately conclude that the function / is analytic at 
(being a solution of analytic ODE), contradicting the assumption that is a point 
of discontinuity of /. 

To prove (fTsT l. we proceed as follows. By Proposition|2]it is easy to see that the term 
Cjj-i ,;• (^)„_i ^1'^"+' is in fact equal to 

in the expression for e^.. Since p„ (^y) ^ 0, we conclude that / {^^^ — f (^/) ^ 0- 
Substituting this into (|9]l, we see that the next term c„_2 j' (^)„-2 ^"^^ equals 

{/'(^;)-/'(^7)}««-2r'"v,(^.), 



and thus /' ^ f ) Proceeding in this manner, we arrive at (ITJi . This 
finishes the proof of Theorem|2] □ 

As a first consequence, we have the real-valued version of the result by Kisunko 

mi. 

Corollary 3. For every n,d > and p ^ consider the family 
^i'ld - |/ e : 2 = E P. W d\ a (D) = d, p„ (x) ^ on [a.b] | . 

Then 

-d-l. 



Since the leading coefficient p„ (x) cannot vanish at more than deg /?„ points, we 
also have the following result. 

Corollary 4. For every n,d > and p ^ consider the family 
=^i5L = |/e^^(2),/^): ^^tpj{x)dj, am=d,degp„<p + 2^. 
Then 

Let us now try to establish what happens in the general case. Let / e ^^(J),/?), 
/ ^ 0. Consider two possibilities. 
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1. The sequence {sk} does not vanish identically. In this case, at least some of its 
initial terms {eq, . . . ,e„(p+2)-i } must be nonzero (Corollary [1}. But this imme- 
diately implies that some of the first n{p + 2) + a—l moments must be nonzero 
as well (otherwise the equality (|6]l cannot hold). 

2. The sequence {e^ } vanishes identically, but Theorem|2]is not applicable. In this 
case the recurrence relation (|6) becomes homogeneous. We rewrite it in the form 

qeik)mk+e = Q, k^Q,l,..., (14) 

e=-n 

where 

= E(-i>'«£+,,;(^+^+;),- (15) 

The leading coefficient qa (k) may have positive integer zeros. Let A (D) denote 
the largest such zero. Then we claim that no more than a+A (D) moments can 
vanish. Indeed, starting with A: = A (S) + 1 we can safely divide the recurrence 
by qa (k) and obtain 

a-i 

m+a = re {k) mk+e, k^A{D) + l, 

t=-n 

where r( (k) are some rational functions with non-vanishing denominators. There- 
fore if the first A (D) + a + 1 moments are zero, then all the rest of the moments 
must vanish, implying vanishing of / itself. 

Thus we have proved the following result. 

Theorem 3. For every 1)^p we have 

(j{D,p) max{«(/:> + 2)- 1, A (S))} + a(2)). 

In Section|6l we demonstrate that in the case of Fuchsian differential operators, the 
number A (£)) has a well-known interpretation. 



5 Moment generating function 

In this section we provide an alternative proof for the result of Kisunko fT2ll con- 
cerning moment generating functions. 

Proposition 5. Let / € The formal power series 

, .def^ Mk 
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is in fact the Laurent series of the Cauchy type integral 

r'f_(t)dt_ 

Proof. Write = 7 (l^) expand into geometric series. Convergence fol- 
lows immediately for °°. □ 

The generalized power sums (Section |3) also have a well-known interpretation as 
the Taylor coefficients of rational functions. The following fact is well-known, and 
so we omit the proof. 

Proposition 6. Let the sequence be of the fonn dlOl ). Then the formal generat- 
ing function 

is a regular at infinity rational function, with poles {^Oj ■ ■ ■ 7 ^p+i }. each with multi- 
plicity at most n. In particular, 

gW = LL (16) 

j=oe=o [z-qj) 

Tlieorem 4. Let f £ ^^(J),/?). Then the Cauchy integral If satisfies in the neigh- 
borhood of °° the inhomogeneous ODE 

^If{z)=Rf{z), (17) 

where Rf{z) is the rational function whose Taylor coefficients at infinity are given 
by the sequence as in Consequently, R f (z) is given by the explicit expression 
il6i , with a( j replaced by en j from (O (Proposition\2}. 

Proof. Consider the asymptotic expansion of the function J) If at infinity 

k=0Z 

By substituting D as in (|2]i and If = Lr=o 3+r we get 

^If = L Pj iz)lf iz) = E L E m, . ^ aij 

j=0 j=0i=0k=Q 
~ 1 " 
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Comparing powers of z we have that Sj^ = /i^ where /i^ denote the left-hand side of 
(|6l). From jU^^ = the conclusion follows. □ 

6 Fuchsian operators 

In this section we employ notions from the classical Fuchsian theory of Unear ODEs 
in the complex domain (we used the reference ifTOl ). 

Assume that the sequence {£<.} vanishes identically. In this case, the Cauchy integral 
If satisfies in the neighborhood of o° the homogeneous ODE 

Definition 5. The operator S) is said to belong to the class if it has at most a 
regular singularity at oo. 

Lemma 2. Let £> G ^H. Then 

1. The numbers (Xj (see Definition^ satisfy 

OCn^Cij, = 0,...,«- 1. (18) 

2. The characteristic exponents ofD at the point °° are the roots of the equation 

qa„ {s-l)=Q, 

where qi (k) is defined by il5i . 
Proof. Dividing the coefficients of T) by p„, we get the operator 

d" + r, (z) d"- ' + ••• + r„ (z) I, o (z) = ^^4^ . 

A necessary and sufficient condition for the point at infinity to be at most a regular 
singularity of this operator is that the function rj (z) is analytic at °o and has a zero 
there of order at least j ( ifTOl Theorem 9.8b]). That is, 

degp„-deg/?„_^- >7. 

But this is equivalent to 

degp„ - n ^ deg/7„_/ - (« - j) 
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To prove the second statement, substitute the formal Frobenius series at infinity 

into D g ~Q. By complete analogy with the calculation in Theorem |4] we get the 
recurrence relation 

f i-^yai.j it+s+i~ i)jbr+.-i =0, f = 0, 1, . . . . 

;■=(), =0 

For r = we find the highest order coefficient in this recurrence to be equal to 
(!-;■ = a = a„) 

n 

£ (-1)-' {s + a + j- l),fl;+a,; =qa{s-l). 
j=0 

The proof is finished. □ 
Together with Theorem|3] this immediately implies the following bound. 

Corollary 5. Let S) G 9^, and let A (2)) denote its largest positive integer character- 
istic exponent at the point °°. Then A (£•) = A (S) — 1, and consequently 

<7{D,p) ^ max{(p + 2)«,A(S))} + t/„-7i- 1. 

Now let us return to Example[T] The following fact is well-known (e.g. IfTOl Section 
9.10]. 

Proposition 7. The Legendre differential operator Dm is of Fuchsian type with sin- 
gularities — 1 , 1 , oo. The characteristic exponents at °° are m + 1 and —m. 

Theorem|2]is clearly not applicable. Using the formula ([7}, it is easy to see that 

Ps(/y)(l)=^2,(/,-^')(-l)=0 

for any / G and therefore the sequence {e^^} in this case is identically zero. 

Consequently, we conclude that 

C7(S„,,0) =OT, 

as expected. 
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